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The possibility of the formation of a condensate of charged spin-0 nuclei 
inside white dwarf cores, studied in arXiv:0806.3692 and arXiv:0904.4267 



is 



C/} Abstract 
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2 pursued further. It has been shown, for cores composed mainly of one element 

(Helium or Carbon), that after condensation phonons become massive and 
, the specific heat drops by about two orders of magnitude. In this note we 

extend that analysis by considering the coexistence of the nuclei of both types 
(Helium and Carbon), whose condensation points are generically different. An 
effective field theory is developed to describe the system when both elements 
are condensed. The spectrum of fluctuations of this two component charged 
condensate possesses a collective massless mode with w oc k . Assuming that 
the fraction of the less abundant element is greater than ~ 1/100, the thermal 
history changes as follows: There is a softer discontinuity in the average specific 
heat after the condensation of first sector, resulting in slower cooling and a 
milder drop in luminosity function. The specific heat remains almost constant 
until the condensation of the second sector, then starts to declines as T 3 / 2 . 
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1 Introduction and Summary 



In a number of previous papers the possibility of forming a charged condensate in 
the core of white dwarfs, composed mainly of a single element (Helium or Carbon), 
and its impact on the cooling process has been studied [HI21E]- It was observed that 
upon condensation the phonons which in standard scenarios of white dwarf evolution 
are the main carriers of internal energy, are eaten by the EM gauge field, become 
massive and freeze out (see also [4]). In this circumstance the relatively small specific 
heat of the degenerate electrons {c v < 10~ 2 ) remains as the only contributor to total 
heat capacity, hence the acceleration of the cooling process. 

In this paper we study the condensation of a white dwarf core containing both 
He and C . This can in particular happen in the following situations: 

1- Carbon density pc is above ~ 10 9 gr/cm 3 , so that the critical temperature for 
condensation, T C C ^J exceeds crystallization temperature T^^Q at about 2 x 10 7 K. 
Pc should also remain below the neutronization threshold ~ 4 x 10 l0 gr/cm 3 . In 
addition we allow the Helium fraction, Xu e = Nn e /N tota i, to range from ~ 1/100 
(when the heat capacity of the impurity starts to dominate that of electron gas) up 
to ~ 1/2. In this range T^ st < T* e . 

2- When the white dwarf is composed mainly of He at densities higher than or 
equal to ~ 10 6 gr/cm 3 , as in low mass Helium white dwarfs (M ^ O.5M ) [TJ [8]. 
We also assume the existence of Carbon impurities greater than 1/100. In this 
case T c c < T^ ryst < T^ e , so one may expect that after the condensation of He 
nuclei the C nuclei form a crystal at T^ ryst . However due to screening effects of 
ife-condensate (studied in [9j[l0]), the C nuclei barely interact with each other and 
fail to crystallize^} 

Therefore it is likely that the He and C nuclei in the core condense as the 
temperature drops below their T c . The condensate can well be described by an 
effective field theory similar to the one developed in [TT1 fl~2] . Nevertheless, after the 
condensation of both He and C, there are two elements involved and in §2 we will 



1 Condensation temperature of a non-interacting Bose gas is given by T\ ni -, — 

27m^ /3 /(£(3/2)) 2/ ' 3 mx, where mx is the mass and nx the density of type X ions. However 
the transition temperature is somewhat higher for interacting particles [5][B], and may be approx- 
imated as the temperature at which the de Broglie wavelengths of the ions of each type start to 
overlap: T^m ~ 47r 2 /3mxd x , where dx — (S/innx) 1 ^ 3 is the average interparticle separation. 
This gives about an order of magnitude higher values than T^ ni -. . 

2 The crystallization temperature is the temperature at which the ratio of coulomb interaction 
of adjacent ions to temperature, T, becomes sufficiently large: T = [(Ze) 2 /4ird}/T crys t ~ 180. 

3 Thus, when the second sector forms a small fraction of the the star, it is reasonable to use 
T*, ni ^ instead of T^u\ as the transition temperature for this sector. To encode this variability we 
define Ax according to 
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lf = A»5_. (i) 
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modify the effective theory to describe it. We will show that the two component 
charged condensate, unlike its one component counterpart, possesses a massless 
collective mode in the spectrum of fluctuations with w oc k 2 . This mode contributes 
to specific heat and changes the thermal history significantly. 

The thermal evolution of the white dwarf, analyzed in §3, would be as follows. 
At first the system can be approximated by an ideal gas with c v = 3/2, and the 
log — log plot of luminosity function (LF) in terms of luminosity has a slope of —5/7. 
Upon the condensation of the first component the corresponding phonons freeze out 
[21 H], thus there is a fall in LF curve whose depth depends on the fraction of the 
second element. The specific heat of the second element's nuclei before condensation 
is approximately temperature independent (= 3/2), resulting in a LF curve with the 
same —5/7 slope. Finally as T c of the second component is approached there is no 
sharp drop in specific heat because of the existence of the massless mode, rather c v 
starts to decrease as T 3//2 , consequently LF curve continuously changes its slope to 
-2/7. 



2 Two Component Charged Condensate 
2.1 Model and Background 

We use Thomas Fermi approximation to treat electrons as a background density 
J° (Ao), and associate to each type of ion a complex scalar field with the appropriate 
charge and mass. In the circumstance of a white dwarf the ions are non-relativistic, 
therefore the most general effective Lagrangian can be written as [TTl IT2"] 
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+ V 2 t^(XW' X-(D' XyX)- l ^^+A^, (2) 

where = — igA^, D' = — ig'A^, and V\(x) and V 2 (x) are polynomial 
functions. 

This system possesses a condensate solution in which the positive charge density 
of $ and X particles is neutralized by the negative charge of electrons, J°. To 
describe the condensate we perform the change of variables: $ = y/m^ae ta , X = 
y /m xX e ^ j an d choose a gauge in which a = 0, to obtain 

C = -7^ + ^(^om^ 2 -i(^) 2 -^V4) ( 3 ) 
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Setting all derivatives except do/3 equal to zero in the equations of motion, the 
background solution must satisfy 



A, = d^ = 0, (4) 

V 1 (0) = V 2 (Q) = 1, (5) 

Jl + J° x + J° e =0, (6) 

with j£ = gm^a 2 , J x = g'm x x 2 ■ (7) 



2.2 Spectrum 

The spectrum of fluctuations around this background contains two transverse and 
one longitudinal modes of a massive gauge field, and a Nambu-Goldstone boson. The 
existence of the latter is easily understood from the symmetries of the problem. The 
original Lagrangian possesses two global U(l) symmetries corresponding to separate 
phase rotations of the two complex scalar fields. They, however, are spontaneously 
broken by the background solution, giving rise to two massless Nambu-Goldstone 
bosons. Out of these two one combination becomes the longitudinal mode of the 
massive gauge field, while the other one remains massless. In the limit m$, mx — > oo 
the dispersion relations simplify to 

to? = k 2 + m\ + m\ . (8) 
2 1 / m\ m\\ k 4 

U L = ~ A \ ^2 + J XT^2 

$ 

1 / m\ m. 



7 ( — 5- + — f —5 o +m\ + ml, (9) 

4 \mj, m x) m i + m 2 



Ug 4 \m 2 x + m%) m\ + ml~ W (10) 



where we have defined 



m l =ga c , m 2 = g'x c . (11) 

Had we started from a relativistic Lagrangian, these dispersion relations would have 
remained unchanged to leading order, except for a k 2 correction to longitudinal 
photon. In addition, two heavy modes corresponding to pair productions of $ and 
X particles would have appeared in the spectrum. 

Note that when the density of one species tends to zero, M approaches the mass 
of that species. This confirms that the particles of this type are efficiently screened 
by the charged condensate of the dominant element, and the massless mode describes 
the low energy excitations of a non-interacting Bose-Einstein condensate, formed by 
these screened particles. 
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3 White Dwarf Cooling 



Since fusion has terminated inside the white dwarfs, they cool down by radiation 
according to [131 IE 



dT 
dt 



L 

'N 



-C(x He A He + x c A c )m u Tj /2 



(12) 



where c v is the average specific heat per nucleus, L the luminosity, N the total 
number of nuclei, Ax and Xx the baryon number and fraction of ion X, m u the 
atomic mass unit, and C a function of atmospheric composition. In this circumstance 
the luminosity function |14j . which is an observable measure of the cooling rate of 
stars, under the assumption of uniform star formation both in space and time takes 
the form 



<p oc 



d\og{L/L & ) 
dt 



oc L 



i(fc+i) 



-i 



(13) 



Here k and n are defined as exponents describing the temperature dependence of 
specific heat and luminosity 



c„ oc T 



L oc T 4/n 



(14) 



n is 8/7, thus, c v which is the single most important variable in white dwarf evolution 
(given the composition of the star), is read easily from the log — log plot of in terms 
of L (LF curve). 

At high temperature it is a good approximation to treat the ions as a classical gas 
with c v = 3/2. This gives the familiar Mestel cooling curve with a = —5/7. As the 
temperature drops below the T c x = max(T c He , T c c ) (see footnote 1) the X nuclei cease 
to contribute to specific heat because of the gap emerging in the phonon dispersion 
relation (as explained in [2]). Nevertheless the specific heat of the other species (Y) 
which is still in the gas state (c„ =3/2) amounts effectively to c v = 3xy/2 per ion. 
Therefore the LF curve drops by log(xy) units. The curve continues with the same 
slope a = —5/7 afterwards. 

As Tj is approached the model of two component charged condensate can be 
applied to the system. Except for the massless mode (10) in the spectrum of fluc- 
tuations, all the other modes have gaps larger than the temperature and are not 
excited significantly. Therefore the specific heat per ion can be written as 



d_ 
df 



n 



d 3 k 



(2tt) 3 e^o/T 



(15) 



where n = N/V is the total number density of ions (including both He and C nuclei), 
and k max is of order of the inverse interparticle separation of the less abundant 
element, above which the effective field theory description breaks down. 
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Figure 1: Schematic behavior of luminosity function. The dashed line corresponds to 
Mestel regime while the solid line exhibits the condensate behavior far from the transition 
point. 



Well below T c y the upper limit of integral can be taken to infinity yielding 

150(5/2) /2MT\ 3/2 M/ , . f2nMT\ 3/2 , . 

and since m# e < M < mc, there is no dramatic change in specific heat right after 
second condensation. Nevertheless it starts to decline as T 3//2 , therefore the LF curve 
after about 



L ( T c )\ 7. fn x \ 7, (m x \ 7, / A 



L{TJ)J 3 \n Y J 2 *\m Y J 2 * \\ 



log ( \ = -log ( Z± ) _ _i og ( — ) + _i og ( ) j (17) 



(see Eq. (jl|) from the first condensation changes slope to a = —2/7. 

We conclude by giving two concrete examples and a schematic diagram of lumi- 
nosity fuction Fig. [l] (note that the typical variations of specific heat at transition 
points may significantly change the behavior near those points). 

As the first example consider a Carbon dominated superdense white dwarf with 
p ~ 10 w gr/cm 3 , xn e — 0.1, and an envelope composition of 



X~0, F-0.9, Z m ~0.1, (1* 



where X, Y and Z m are H fraction, He fraction and metallicity, respectively. Carbon 
nuclei first reach their condensation temperature at L ~ 3L Q when LF curve drops 
for log(l/x# e ) = 1 unit, and the second condensation takes place later at L ~ 
3 x l(T 4 L m . 
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Next consider a Helium dominated white dwarf with p ~ 10 6 gr/cm 3 , Xc ~ 0.1, 
and atmospheric composition of 

X-0.99, F~0, Z m ~2xl(r 4 . (19) 

In this case He condenses first at L ~ 1.5 x 1CT 4 L with again a unit drop in LF 
curve, and then C nuclei condense when L ~ 6 x 10 _12 L Q . 
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